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Abstract. We study the spectral statistics of interacting spinless fermions in a two-dimensional disordered 
lattice. Within a full quantum treatment for small few-particle-systems, we compute the low-energy many- 
body states numerically. While at weak disorder the interactions reduce spectral correlations and lead to 
localization, for the case of strong disorder we find that a moderate Coulomb interaction has a delocalizing 
effect. In addition, we observe a non-universal structure in the level-spacing distribution which we attribute 
to a mechanism reinforcing spectral correlations taking place in small systems at strong disorder. 

PACS. 71.27.+a Strongly correlated electron systems; heavy fermions - 73.20.Jc Derealization processes 
- 72.15.Rn Localization effects (Anderson or weak localization) 



1 Introduction 

The ongoing miniaturization of electronic devices and 
the peculiar physics associated motivate the investiga- 
tion of systems with reduced dimensionality. In such sys- 
tems, the effect of Coulomb interactions is expected to be 
strong, and prominent experimental observations of the 
last decade are thought to be due to correlation effects. 

One such experimental result concerns the magnitude 
of the persistent currents in disordered mesoscopic rings 
[T|P| , which are much greater than theoretical predictions 
from approaches neglecting correlations |3]- The interac- 
tion effect seems to be beyond the perturbative regime, 
but a full treatment of the realistic situation is out of 
reach (for a review see 0]). Nevertheless, analytical jS] 
and numerical |^||7) calculations in ID and in 2D have 
shown that a repulsive interaction can enhance the persis- 
tent currents. 

As another important example, a metallic behaviour 
' has been observed [H] (see |H] for a review) in 2D elec- 
tron gases at low electronic density, where the ratio r s of 
Coulomb to kinetic energy is large (> 10). This cannot 
be explained by the standard scaling theory of localiza- 
tion which, neglecting electronic correlations, predicts an 
insulating behaviour in 2D for any finite disorder strength 
|l(Jj . Since interactions become important when the elec- 
tron density is low (large r s ), they have been suggested 
to be responsible for the observed metallic behaviour [7|. 
However, a perturbative introduction of the interaction 
leads to a reinforcement of the electron localization |11| . 
This points towards the necessity of treating the interac- 
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tions beyond perturbation theory, and several approaches 
have been proposed in this direction. 

Renormalization group techniques indicate that the in- 
teraction permits a metal-insulator transition in a weakly 
disordered two-dimensional electron gas |12II13| . However, 
despite some well reproduced properties (such as the effect 
of a magnetic field), a description of the transition is not 
possible and the metallic phase is based on the assumption 
of a Fermi liquid. 

Alternatively, field theory based on the fact that the 
metallic phase is not a Fermi liquid (as argued in |14j i 
has shown that a perfect metal can be stable in 2D if the 
interaction is strong enough 15 . 

Furthermore, numerical calculation of the current- 
current correlation function using a quantum Monte Carlo 
approach has shown that interactions change the be- 
haviour of the conductivity at low temperature from an 
insulating to a metallic one |16| , as observed in the exper- 
iments. 

All these approaches contain approximations. On the 
contrary, numerical studies of model systems allow for ex- 
act results, although they are limited to simplified models 
and small system sizes. Therefore, such studies are able to 
provide a complementary view of the mechanisms implied 
in the interplay between disorder and interaction. 

In the present work, we numerically investigate inter- 
acting spinless fermions (spin-polarized electrons) in small 
two-dimensional lattices with disorder, performing a direct 
diagonalization of the Hamiltonian. 

The limitation in size makes it difficult to vary directly 
the electronic density as in experiment. To change the in- 
teraction parameter r s , we vary instead the interaction 
strength U while keeping the system size L and the par- 
ticle number N constant. 
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Table 1. Summary of the behaviour of the distribution P(s) 
for different regimes of interaction U and disorder in small 
two-dimensional systems. PWMS stands for Pinned Wigner 
Molecule Statistics (see text). The cross-over between the weak 
and strong interaction regimes represent the main result of the 
paper. 



U 



Weak disorder Strong disorder 





oo 

cross-over 



Wigner-Dyson 

PWMS 

monotonic 



Poisson 
PWMS 
non- monotonic 



Within this approach, we explore the interaction ef- 
fects on the ground state structure and the probability 
density P(s) of the normalized level spacing s — A/(A), 
where A = E\ — Eq is the energy spacing between the 
many-body ground-state and the first excited state. We 
denote by (...) the average over the ensemble of disorder 
configurations. 

In the non-interacting case, A is equal to the one-body 
level spacing at the Fermi energy. The statistics of these 
one-body level spacings P(s) has been extensively studied 
and found to be an indicator of the metal-insulator Ander- 
son transition occurring in 3d as a function of the disorder 
strength [17I18| . In the diffusive regime, P(s) corresponds 
to Wigner-Dyson statistics, while the Anderson insulating 
regime is characterized by Poisson statistics (table ^1 . 

In the opposite limit, at U = oo, the structure of 
the ground state is imposed by the Coulomb repulsion, 
which leads to a Wigner crystal pinned by the disor- 
der. We show in section Q] that the resulting distribution 
Pu=oa (s) (which we call "Pinned- Wigner-Molecule Statis- 
tics" , PWMS) is non-universal for finite size systems. 

For weak disorder, we found the dependence of P(s) 
on U to be consistent with previous studies: P(s) crosses 
over smoothly from Wigner-Dyson statistics to its infinite 
interaction limit 19 20.21 . 

Within an approximative method (Configuration In- 
teraction method, starting from Hartree-Fock orbitals), 
the strong disorder case has been studied by Benenti et al. 
for larger systems 22 . They obtain an interaction-induced 
transition of P(s) from Poisson to Wigner-Dyson. Our 
findings of an increase of spectral correlations by mod- 
erate interaction show that this behaviour persists when 
electronic correlations are taken fully into account. 

The non-trivial behaviour of the spectral statistics for 
strong disorder is discussed in section [3 after the presen- 
tation of the model we study (section |2J). In section 0J we 
present results for the inverse participation ratio of the 
ground state in the many-body on-site basis, and we show 
that the increase of spectral correlations is related to a de- 
localizing effect. Our conclusions are discussed in section 
03 and a small size effect modifying the shape of P(s) is 
presented in the Appendix. 



2 Disorder and Coulomb Interaction 

We consider TV spinless fermions on a two-dimensional L 
by L lattice (in the following, we concentrate on N = 4 
and L = 6). We note M = L 2 the number of sites. 

The Hamiltonian of the system is H = Ha+Hjj, where 
Ha is the standard Anderson Hamiltonian 



H A = -*£(< 



M 



Cj c 



>+E' 



(1) 
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with cf (cj) creating (destroying) an electron on the site i. 
The first term of Ha allows for hopping between nearest 
neighbours < i,j > on the lattice. We take t = 1, rep- 
resenting then the energy scale. Periodic boundary con- 
ditions are used, leading to a toroidal topology. The sec- 
ond term of Ha models a disorder potential. The Vi are 
independent random variables uniformly distributed in 
[-W/2; W/2], and W is the disorder strength. 

The interaction term is chosen to be of the Coulomb 
form: 

H - U V^ C J C% Ci Cj (o\ 

fft/ -iL z: ' (2) 



where dij is the smallest distance on the torus between 
the sites i and j and U is the interaction strength. 

This model allows us to study qualitatively the effect 
of interactions in disordered systems. For the exact diag- 
onalization of the Hamiltonian H, we have used a routine 
developed by Simon and Wu (22| based on the Lanczos 
algorithm |24j . 



3 Interaction induced many-body 
level-repulsion 

In this section, we study in detail the probability distri- 
bution P(s) of the first excitation energy whose behavior 
in different regimes is sketched in tabled 

In the absence of interactions, Pu=o(s) corresponds to 
the Wigner-Dyson distribution Pwd(s) = ^sexp(— js 2 ) 
for a diffusive (metallic) system. For a strongly disordered 
system (Anderson insulator), its limiting behavior for L — * 
oo is the Poisson distribution Pp(s) = exp(— s) 1 . 

In the strong interaction limit (U — ► oo), the elec- 
trostatic energy dominates and the energetically lowest 
many-body states correspond to periodic distributions of 
the electrons (Wigner crystal) which are pinned even by 
a small amount of disorder. All Wigner crystals have the 
same electronic structure, thus the same Coulomb energy, 
but differ in their location on the lattice and therefore in 
their disorder energy. 

In the limit of low densities, the ground state struc- 
ture approaches the usual hexagonal Wigner crystal, the 



1 In finite systems, Pu=o(s) is distinct from Poisson even for 
W — > oo. For the system size we consider (M = 36), however, 
the difference is rather small. 
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number of possible positions on the lattice is large, and 
Pu=oo(s) tends to Poisson. For higher densities we need to 
take into account commcnsurability effects. In the studied 
case, N and L are such that the nine energetically low- 
est U = oo many-body states are square-shaped Wigner 
crystals (which we will refer to as "Wigner molecules") 
differing only by their location on the lattice. These nine 
Wigner molecules have different disorder energy, and the 
first excitation energy is 



= E X - 


N 

fc=i 


N 

-£ 

fc=i 



J 3k i 



(3) 



where the jk describe the sites occupied by a particle 
in the energetically lowest Wigner Molecule (the ground 
state) and the ik the occupied sites in the second Wigner 
molecule (the first excited state). Since we only have 
nine configurations, Pu=<x>{s) is not exactly Poisson. In- 
stead, it is a distribution which is intermediate between 
semi-Gauss 2 and Poisson, which we call pinned- Wigner- 
molecule statistics (PWMS). 

Figure QJi shows P(s) with W = 5, for three values 
of U. With this disorder strength, the single-particle lo- 
calization length is larger than the system size, so that 
the motion of non-interacting particles through the sam- 
ple is diffusive and P(s) follows Wigner-Dyson statistics 
at U = 0. 

In this case, we find that the level repulsion decreases 
with increasing interaction strength. In this sense, this 
result is consistent with those of 19,20,21 . The difference 
is that the U — oo limit of P(s) is in fact the PWMS and 
not Poisson for the investigated system sizes. 

This behaviour is as expected: strong interactions drive 
the system out of the diffusive regime by setting up a 
Wigner Crystal which is pinned by disorder. Therefore 
spectral correlations disappear with increasing interaction 
strength. 

Figure ^ shows our results for the case of stronger 
disorder (W = 20). This disorder is not strong enough 
to have a Poissonian P(s) in the non-interacting case, 
and hence Pu=o(s) is intermediate between Poisson and 
Wigner-Dyson. Nevertheless, a single-particle state at the 
Fermi energy is localized over about three sites only, mean- 
ing that the system is far from the diffusive regime for 
U = Q. 

The main feature is that at this disorder strength, the 
evolution of P(s) with the interaction strength is non- 
monotonic. Obviously, for stronger interactions the spec- 
tral correlations eventually decrease and P(s) approaches 
its infinite interaction limit (PWMS). For a moderate U, 
however, the spectral correlations are increased relative to 



2 At half filling there are only two different Wigner crystal 
configurations. Their energy difference is given by the differ- 
ence between two sums of independent random numbers. Ac- 
cording to the central limit theorem, its distribution in the 
limit of N — > oo is Gaussian. Since only positive values for s 
are considered, the resulting distribution P(s) is called semi- 
Gauss. 
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Fig. 1. Distribution of the lowest many-body excitation for 
different interaction strengths U at W — 5 (a) and W = 20 
(b), for 4 particles on a 6 by 6 lattice. Each curve is obtained 
from data for 9000 disorder configurations. 



the non-interacting case (i.e. P(s) is more Wigner-Dyson- 
likc). This is the main result of the present paper. The 
increase of spectral correlations could be the precursor of 
the transition towards universal correlations found in |22j . 

The non-monotonic behavior becomes less significant 
as W increases further. It is for this reason that we have 
presented our results for W = 20. Note that, with inter- 
action, a peak appears in the P(s) curves. This peak is 
the manifestation of a mechanism enhancing spectral cor- 
relations in a non-universal way for small system size. We 
discuss this mechanism in detail in the appendix. 

In order to characterize quantitatively the non- 
monotonic behaviour exhibited by P(s), we consider the 
evolution of the variance of s with U . We show in figure|21 
that this evolution in the cases of W = 5 and W — 20 is 
very different. At W — 5 there is a monotonous evolution 
towards smaller correlations, while at W = 20 the variance 
exhibits a minimum as a function of U (corresponding to 
an increase of the spectral correlations) . 

With a short range interaction, we have found sim- 
ilar results for P(s) at half filling. However, away from 
half filling, configurations not affected by interactions are 
connected by hopping matrix elements. Therefore, short 
range interactions cannot suppress the mobility of the elec- 
trons. As a consequence, the variance (or another param- 
eter characterizing the distribution P(s)) does not reach 
its PWMS value at strong U. 
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Fig. 2. The variance of s as a function of U at W = 5 (filled 
circles) and W = 20 (open squares), for 4 particles on a 6 by 
6 lattice. Each point is computed from 3000 disorder configu- 
rations. The statistical error is smaller than the symbol size. 
The W = 20 infinite U limit PWMS is approached from above 
at values of U which lie outside the scale of the figure. For 
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4 Delocalization in the many-body on-site 
basis 

The fact that P(s) approaches the Wigner-Dyson distri- 
bution in the presence of a moderate interaction could 
be interpreted as a signature of a delocalization of the 
electrons. Whereas the link between spectral correlations 
and localization is clear for one-particle level statistics 
(U = 0), it is less obvious in the interacting regime. To 
clarify this interpretation, we have studied the localization 
of the ground-state 



is) = 5>»i»> 



(4) 



in the many-body on-site basis {|n)}, via its inverse par- 
ticipation ratio 

(5) 



n-* = Y,\*n 



Contrary to P(s), the inverse participation ratio TiT 1 
depends on the choice of the basis in which it is calculated. 
To be allowed to interpret the inverse participation ratio 
as a measurement of the electron localization, we have 
chosen for {\n}} the Slater determinants 



C ii C i 2 C i 3 C U 



10} 



(6) 



which correspond to the four particles being localized on 
the lattice sites i\, 12, *3 and 14 and where |0) is the empty 
lattice state. 

If \g) is given by one of these basis states \m), then 
&n = 5n,m, and "Rr 1 = 1. If, on the other hand, \g) is 
a superposition of many elements of {|rc)}, then IZ^ 1 is 
very small compared to unity. Since the basis is built with 
completely localized electrons, "Rr 1 can be interpreted as a 
measurement of the localization of electrons in the many- 
body state \g). 



Fig. 3. Evolution of 1/ (72T 1 ) as a function of U for W = 5 
(circles, in the inset), W = 20 (squares), W = 23 (diamonds), 
W — 25 (triangles) and W = 30 (crosses) . The statistical error 
is smaller than the symbol size. 



Figure |3 shows numerical results for the inverse of the 
average of Hr 1 . In the case of weak disorder (W = 5, 
in the inset), this quantity is monotonically decreasing 
with U, as expected. The U = diffusive situation with 
delocalized particles is perturbed by moderate interactions 
which increase the scattering, reducing both the mobility 
and 1/ (TZ^ 1 '). In the regime of strong interactions, the 
electrons form a Wigner crystal to minimize electrostatic 
energy. Since this corresponds to one particular state of 
the chosen basis, Tlr 1 decreases to one. 

In the regime of strong disorder, the behaviour is very 
different. At weak interaction, 1/ \JZ~] increases with 
U, which means that interaction has a delocalizing effect. 
This can be understood within the following scenario. 

At U = 0, we have localized one-body wave functions, 
the disorder being strong enough to dictate the electronic 
configuration. For finite U, the effect of the interaction de- 
pends on this particular sample-dependent configuration. 
In some samples, this electronic structure is close to the 
one adapted to interaction (Wigner molecule), therefore 
the interaction strengthens the localization. On the other 
hand, in most samples, the disorder-adapted electronic 
configuration is rather different from the Wigner molecule 
structure. Therefore, increasing interaction strength in- 
duces charge reorganizations (as proposed in 25 in the 
context of one-dimensional rings) at particular sample- 
dependent values of the interaction strength U c . When 
U ~ U c , the competition between interaction and disor- 
der leads to a ground state which is a superposition of 
a state adapted to disorder and another one adapted to 
interaction. This results in a pronounced delocalization. 

Since U c is strongly sample-dependent, these dereal- 
izations are smoothed by the disorder average. Further- 
more, we expect that as W increases, the mean value of 
U c also increases and its distribution spreads, consistent 
with the behaviour observed in figure 

To illustrate these considerations, we can define for 
each sample the increase 



5 = K- 1 (U)-K Q 



(7) 
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Fig. 4. Distribution of S at W = 20 for U = 2.5 (filled squares) 
and U — 9 (open squares) calculated from 3000 disorder con- 
figurations. The averages are indicated below the curves. The 
statistical error for these averages is smaller than the symbol 
size. 



of TlT 1 with respect to the non- interacting value TZ^ 1 . Fig- 
ure 0] depicts the distribution of this quantity at W = 20, 
for two different values of U. We can see that the most 
probable value of S is close to zero. Nevertheless, for 
U = 2.5 the negative tail of the distribution, correspond- 
ing to strong delocalizations, is more pronounced than the 
positive one, which corresponds to strong localizations, re- 
sulting in a negative average of <5. This is not the case when 
the interaction is strong (U = 9). We can conclude that 
for W — 20, U c is more probable to be close to 2.5 than 9. 



5 Conclusion 



We have studied the spectral statistics of interacting spin- 
less fermions in a two-dimensional disordered system, us- 
ing exact diagonalization. 

We have found that correlations in the statistics of the 
level spacing between the many-body ground state and the 
first excited state are increased by a moderate interaction 
when the disorder is strong enough to localize the one- 
particle wave functions. The interpretation of this effect 
as a delocalization of the electrons has been supported by 
the study of the inverse participation ratio for the ground 
state. 

The delocalization effect can be understood as a conse- 
quence of a competition between disorder and interaction 
for the structure of the ground-state, taking place at some 
sample- dependent value of the interaction strength. This 
competition results in charge reorganization caused by in- 
creasing interaction from a configuration minimizing the 
disorder energy to a configuration more adapted to the 
interaction. 



We acknowledge very useful discussions with and comments by 
R. A. Jalabert, P. E. Falloon, G.-L. Ingold, and J.-L. Pichard. 



A Appendix: Non-universal level statistics in 
small systems 

The peak appearing in the P(s) curve of figure \I]p for 
U = 20 is the manifestation of a mechanism enhancing 
spectral correlations in a non-universal way. Since its ma- 
jor ingredient is a competition between disorder and in- 
teraction, this mechanism could be the precursor of what 
happens in bigger systems, where other mechanisms could 
take place at higher order, resulting in enhanced spectral 
correlations even in the thermodynamic limit. 

The mechanism we want to describe takes place in 
very small systems when the mean level spacing (A) is 
not much smaller than the hopping t. In order to explain 
this mechanism, we start by considering the simpler non- 
interacting situation, before treating the interacting case. 



A.l Non-interacting case 

Without the hopping t, the eigenstates of the system are 
Slater determinants of particles localized on single sites, 
and the distribution of their energies is Poisson. In the 
limit of strong disorder, W ^> t, the coupling t of neigh- 
bouring sites is typically much smaller than their differ- 
ence in on-site energy. Therefore, particles remain local- 
ized on single sites, except in samples for which the energy 
of the highest occupied site i is almost degenerate with 
that of the lowest unoccupied site j, and with these two 
sites being nearest neighbours. 

In those samples, the hopping couples directly two al- 
most degenerate levels, resulting in a delocalization of one 
electron over these two sites (as pointed out in |2l)]1. 

If the mean spacing (A) is bigger than t (small size or 
very big W), the two coupled levels can be considered as 
a two-level system. Therefore the level spacing is 



A 



\J{Vj - Vi 



At 2 . 



(8) 



In the case where the sites i and j are not nearest neigh- 
bours A — Vj — Vi ~ 0. As a consequence, the special 
samples are responsible for the appearance of a dip and a 
peak in P(s) at s = and s = 2t/ (A), respectively. 

On the other hand, if the mean spacing (A) is smaller 
than t (greater size and not too big W), typically more 
than two levels are coupled and eventually universal 
random-matrix-theory-like correlations can arise. Even 
though at intermediate sizes the spectral correlations are 
still greater in the special samples, their weight in the en- 
semble rapidly decreases with L. In the thermodynamic 
limit, the anomaly therefore disappears. 



A. 2 Interacting case 

In the absence of hopping (t — 0), given a quite strong 
disorder W, it is always possible to find values of U such 
that the two energetically lowest many-body states of a 
given sample are almost degenerate. For example, one can 
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minimize the interaction energy (with a Wigner crystal) 
and the other one can be slightly different, increasing the 
interaction energy while reducing disorder energy. 

Now, the first two many-body states can in certain 
samples be connected by only a single hop of one particle. 
The probability for such a situation is quite large since 
one many-body state is coupled to many others by the 
hopping of one of the particles. 

In those samples, the introduction of the hopping cou- 
ples directly two almost degenerate levels. If the values of 
U and W are strong enough, the particles remain localized 
on individual sites except the one implied in the connec- 
tion of the two states. This results in a delocalization of 
the ground state in the on-site basis (1Z = 2). 

For (A) > t (small size or very big W) , the two coupled 
levels can again be considered as a two-level system, and 
the non-universal correlations appear as explained previ- 
ously. 

For (A) < t (greater size and not too big W), typically 
more levels are coupled and the two-level system approxi- 
mation breaks down. As for the non-interacting case, this 
mechanism disappears in the thermodynamic limit. How- 
ever, other mechanisms involving more hoppings can take 
place, and the competition between disorder and interac- 
tion might still induce delocalization and eventually stem 
the spectral correlations seen in j22j . 



15. S. Chakravarty, L. Yin, and E. Abrahams, Phys. Rev. B 
58, R559 (1998). 

16. P. J. H. Denteneer, R. T. Scalettar, and N. Trivedi, Phys. 
Rev. Lett. 83, 4610 (1999). 

17. B. L. Al'tshuler and B. I. Shklovskii, Zh. Eksp. Teor. Fiz. 
91, 220 (1986), [Sov. Phys. JETP 64, 127 (1986)]. 

18. B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, 
and H. B. Shore, Phys. Rev. B 47, 11487 (1993). 

19. G. Benenti, X. Waintal, and J.-L. Pichard, Europhys. Lett. 
51, 89 (2000). 

20. R. Berkovits, Y. Gefen, I. V. Lerner, and B. L. Altshuler, 
Phys. Rev. B 68, 085314 (2003). 

21. G. Katomeris, F. Selva, and J.-L. Pichard, Eur. Phys. J. 
B 31, 401 (2003). 

22. G. Benenti, X. Waintal, J.-L. Pichard, and D. L. Shep- 
elyansky, Eur. Phys. J. B 17, 515 (2000). 

23. H. D. Simon and K. J. Wu, TRLan software pack- 
age (http://www.nersc.gov/research/SIMON/trlan.html, 
1999). 

24. J. K. Cullum and R. A. Willoughby, Lanczos algo- 
rithms for large symmetric eigenvalue computations, vol. 1 
(Birkhauser, Boston, 1985). 

25. D. Weinmann, P. Schmitteckert, R. A. Jalabert, and J.-L. 
Pichard, Eur. Phys. J. B 19, 139 (2001). 

26. A. Wobst, G.-L. Ingold, P. Hanggi, and D. Weinmann, 
Phys. Rev. B 68, 085103 (2003). 



References 



1. L. P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, Phys. 
Rev. Lett. 64, 2074 (1990). 

2. V. Chandrasekhar, R. A. Webb, M. J. Brady, M. B. 
Ketchen, W. J. Gallagher, and A. Kleinsasser, Phys. Rev. 
Lett. 67, 3578 (1991). 

3. A. Schmid, Phys. Rev. Lett. 66, 80 (1991); F. von Oppen, 
E. K. Riedel, ibid., p. 84; B. L. Altshuler, Y. Gefen, Y. 
Imry, ibid., p. 88. 

4. U. Eckern and P. Schwab, Advances in Physics 44, 387 
(1995). 

5. T. Giamarchi and B. S. Shastry, Phys. Rev. B 51, 10915 
(1995). 

6. E. Gambetti-Cesare, D. Weinmann, R. A. Jalabert, and 
P. Brune, Europhys. Lett. 60, 120 (2002). 

7. G. Benenti, X. Waintal, and J.-L. Pichard, Phys. Rev. Lett. 
83, 1826 (1999). 

8. S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, 
V. M. Pudalov, and M. D'lorio, Phys. Rev. B 50, 8039 
(1994). 

9. S. V. Kravchenko and M. P. Sarachik, Rep. Prog. Phys. 
67, 1 (2004). 

10. E. Abrahams, P. W. Anderson, D. C. Licciardello, and 
T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 

11. B. L. Altshuler, A. G. Aronov, and P. A. Lee, Phys. Rev. 
Lett. 44, 1288 (1980). 

12. C. Castellani, C. Di Castro, and P. A. Lee, Phys. Rev. B 
57, R9381 (1998). 

13. Q. Si and C. M. Varma, Phys. Rev. Lett. 81, 4951 (1998). 

14. V. Dobrosavljevic, E. Abrahams, E. Miranda, and 
S. Chakravarty, Phys. Rev. Lett. 79, 455 (1997). 



